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Abstract 

We construct fractionally integrated continuous time GARCH models, which capture the observed 
long range dependence of squared volatilities in high-frequency data. We discuss the Molchan-Golosov 
kernel and the Mandelbrot-van-Ness kernel with respect to their integrability, and the resulting fractional 
processes with respect to their increments, as they should be stationary and positively correlated. Since 
this poses problems we resort to moderately long memory problems by choosing a fractional parameter 
d 6 (—0.5,0) and remove the singularities of the kernel to obtain non-pathological sample path properties. 
The new fractional COGARCH process has certain positive features like stationarity and algebraically 
decreasing covariance function. 


1 Introduction 

The goal of this paper is to construct fractionally integrated continuous time GARCH (COGARCH) models 
in order to capture the long range dependence behaviour of squared volatilities as observed in high- 
frequency data in finance, but also in turbulence data and others. All stochastic quantities of our paper are 
defined on a probability space (O, T, P). 

COGARCH models are specified by two equations, the mean and the variance equation. The single 
source of variation is a Levy process L with characteristic triplet ( jt, a 2 ,Vi). We refer to fl20| for back¬ 
ground on Levy processes. Then the COGARCH(l,l) price process ( Gt)t>o is given by 

dGt — Ct-dLt, t > 0, 


where the squared volatility a} satisfies the SDE 

da} = —fii (a}_ — ao)dt + ot\o}_d[L, L]^°\ t > 0, (1.1) 

with parameters <x.q,oc\,^>\ > 0. When [L,L]t denotes the quadratic variation process of L, then 

[L, L\t = a 2 t + A L s ) 2 = a 2 t+[L,L} { t D) , t > 0, 

0<s<t 


and [L, L]j denotes the discrete part of the quadratic variation process of L. 
The solution of m is given by 


= e~ Xt 
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■«O01 [ < 
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S ds 


t > 0, 


( 1 . 2 ) 
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with auxiliary process 


Xf = |6if-X>g(l + «i(AL s ) 2 ), f >0. (1.3) 

S<t 

For more background on COGARCH processes cf. HI and for a recent review |9]- The auxiliary process 
X is a spectrally negative Levy process driven by the subordinator [L, L]* D k In order to obtain any long 
range dependence model for the volatility this driving subordinator has to be modified. 

Fractional Levy processes (fLp) generalize fractional Brownian motion (fBm) in a natural way. It is 
well-known (e.g. j5j or fl9l ) that fBm can be defined as a convolution of Brownian motion with a Volterra- 
type kernel. Two such kernels with fractional parameter d G (—0.5,0.5) are the Mandelbrot-van-Ness 
kernel, which leads to fBm on R and the Molchan-Golosov kernel, which results in fBm on R+ = [0, oo). 
Such Gaussian processes have continuous sample paths, stationary increments, and they are self-similar. 
Moreover, they can model long range dependence for d > 0. 

For fractional parameter d G (—0.5,0.5) the Molchan-Golosov (MG) kernel is defined for t > 0 as 

G (f,s) = c d {t - s) rf 2 Fi (~d,d,d + 1, , S G [0,f], (1.4) 

where 2 -Fl is Gauss' hypergeometric function, and the Mandelbrot-van-Ness (MvN) kernel is defined for 
t G If as 

ff vN (t,s) = c d ((f - s)i - (s)i) , set, (1.5) 

with constants c d = ( rfi+d)ri '{i-2d ) ) ~ anc ^ Td = r(d+i) ' res P ect i ve ly- Replacing the driving Brownian 
motion by a non-Gaussian Levy process leads to fractional models with different finite-dimensional distri¬ 
butions (cf. El)- 

Both representations are useful depending on the envisaged application. For instance, the Molchan- 
Golosov fLp has no infinite history, whereas the Mandelbrot-van-Ness fLp has. More properties have been 

shown e.g. in to, m m, and El- 

According to 113 the existence of a fractional subordinator integral with respect to a Levy process L re¬ 
quires the following properties: L has finite second moment and the kernel / G L 1 flL 2 . Unfortunately, the 
MG kernel leads to a fractional subordinator with non-stationary increments and the MvN kernel belongs 
to L 1 D L 2 only for negative fractional parameter. With the goal to obtain non-pathological sample paths 
we define a modified MvN kernel. Using this kernel function we obtain a fractional subordinator, which 
has a continuous modification and stationary increments. The autocovariance function of these increments 
decreases with an algebraic rate. This allows us to define a fractionally integrated COGARCH(l,l) volatil¬ 
ity process driven by the subordinator [L,L] I ' D ' 1 . We also find a stationary version of this new volatility 
model, which results in a price process with stationary increments. 

Our paper is organised as follows. We will discuss fractional subordinator models based on the MG-fLp 
and the MvN-fLp for d G (—0.5,0.5) in Section [2] In Section[3]we present a new modified MvN kernel and 
investigate its properties. We also define the resulting long memory subordinator and present the sample 
path properties, its cumulant generating function, and the properties of its increments. Taking as subor¬ 
dinator the quadratic variation process of a finite variance Levy process allows us to define a fractional 
COGARCH(l,l) process in Section [4] Our main result in this section is the existence of a stationary version 
of the variance process, which implies that the price process has stationary increments. We conclude our 
paper with the straightforward extension of our approach to COGARCHfp, q) processes. 


2 Fractional subordinators - naive approaches 

The following result will clarify the existence of the fractional subordinator driven integral. 
Proposition 2.1. Consider the kernels in m and O- Then the following holds. 
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(i) If de (-0.5,0.5) then / d MG (f,-) £ L X (R) n L 2 {R)for all t > 0. 

(ii) Let t £ R. If d £ (0,0.5), then f^ vN (t, ■) £ L 2 (R) but f™ vN (t, •) g L 1 (R). 

I/, however, d £ (—0.5,0) f/zen f^ vN {t, •) £ L 1 (R) n L 2 (R). 

Proof. For part (i) see 13 Remark 3.3]. For the proof of (ii) see 13 Proposition 2]. □ 

Consequently, (cf. |17 Theorem 3.3]) we can define fractional subordinators based on the MG kernel 
for all d £ (—0.5,0.5) as 


Sf = fff G (t,u)dS u , f>0. (2.1) 

Proposition 2.2. Let S — (St)t>o be a subordinator satisfying E(S 2 ) < oo. Then Sf as defined in d2.1b exists for 
all t > 0 as limit in probability or in the L 2 (d)-sense. 

Fractional subordinator models have also been suggested in |[2j for a time change process on [0,oo). 
Their Example 1 is based on d > 0. In this case / f ') /,G (-, s) is increasing for each s > 0 and since they require 
that S is a strictly increasing subordinator, S d is a.s. strictly increasing. Hence the integral d2.Hl exists as an 
improper Riemann integral. Further, S d is adapted to the filtration generated by S such that 1161 Section II, 
Theorem 7] implies that S d is a semimartingale with characteristics (S rf , 0,0). 

Since we would like to model the long range dependence effect seen in empirical volatilities, we are 
interested in the case d > 0. However, a simple calculation shows that for every d £ (—0.5,0.5) \ 0 the 
expectation E(Sf) = C,/f 2rf ] , which implies immediately that S' 1 cannot have stationary increments for 
d 0. This makes S d as driving process of any SDE awkward, since the calculation and estimation of the 
model parameters would be extremely hard, if at all possible. 

Consequently, we turn to the fractional subordinator based on the MvN kernel, which requires a two- 
sided subordinator. As usual, for any Levy process L, we define a two-sided Levy process 

Lt = -L ( lll {i<0} + L?h m} ( 2 . 2 ) 


on the whole of R, where L ! | ' and L 121 are two independent and identically distributed copies of L. Then 
for a two-sided subordinator S we define for d £ (—0.5,0) (where the integral exists) 


Sf= [ f^ vN (t,u)dS u , f£R. 

J IR 


(2.3) 


However, having the same covariance structure as fBm, S d has negatively correlated increments for 
d £ (—0.5,0). Further, due to the singularities at s = 0 and s — t the kernel f^ vN (-,s) is discontinuous 
and unbounded for all s £ R, such that S l1 has discontinuous and unbounded sample paths with positive 
probability (cf. 1211 Remark 3.3] and [18l Theorem 4]). 


3 Modified MvN subordinators 

There is still a possibility to construct some long-range dependent fractional subordinator with stationary 
increments following an idea of Brockwell and Marquardt fl3l . Section 8. In order to define fractional 
CARMA subordinator models, they suggest a convenient family of moderately long memory fractional 
Levy processs generated from the MvN kernel f^ /lvN in 111. 5b . We will adapt their approach in principle. 
First we consider the MvN-kernel / ( ^ tcN (f, •) for t £ R and restrict d to be negative, that is 

-(-<). S6R. (3.1) 

From Proposition 12.11 we know that f^ vN (t, •) £ L 2 (R) fl L 2 (R). 
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Figure 1: Simulated MvN-fLp (dashed line) with fractional integration parameter d — 0.25 (left) and 
d — —0.25 (right) driven by a compound Poisson process (solid line) with rate 0.4 and standard normally 
distributed jump size. 


A fractional Levy process generated from the kernel J3.H l will still have the same covariance structure 
as fractional Brownian motion. Hence the increments of such a process will be negatively correlated, see 
e.g. nu Section 5.1]. Further, due to the singularity at s = f the function , s) is discontinuous for 

all s £ R, such that S cl has discontinuous and unbounded sample paths with positive probability (cf. Iil8l 
Theorem 4]). To overcome these drawbacks, we will bound the kernel at its singularities. 

Observe that the MvN-kernel f^ vN (t, ■) is up to a constant given by s H > gd{t ~ s ) ~ gd(~ s ) where the 
function gj is defined by gd(x) := xi. This suggests to bound gd at the singularities s = 0 and s = f by 
incorporating a shift a > 0, leading to 

gaA x ) := ( a + x +) d ' xeR. (3.2) 

In the following we give the definition of the resulting modification of the Mandelbrot-van-Ness kernel. 
Definition 3.1. Let d < 0 and a > 0. For each t 6 R the non-normalized modified MvN-kernel is given by 

/«,d( f ' s ) = ga,d(~ s ) ~ gaA* - s) 

- (a + (-s) + ) d - (fl + (f-s) + ) rf , seR. 

Remark 3.2. (a) The normalizing constant of the MG kernel is needed for the integral representation offBm, and 

for the MvN kernel the constant is a natural consequence of fractional integration or differentiation. As this is 
not obvious for our new kernel, we work with the unnormalized version. 

(b) Observe that besides substituting g^ by g a ^, we changed the signs. This way we ensure f a ,d(t, •) to be non¬ 
negative for t > 0. 

Proposition 3.3. For a > 0 and d < 0 consider the modified MvN-kernel f a $ as in J3.3I I. Then the following holds 
for all t G R. 

(i) f a/ d{t, ■) is continuous, 

(ii) |/ fl/ d(f,s)| < a d for all s £ R, 

(iii) |/ fl( d(f,s)| ~ |s| d_1 as s -> -oo. 

(iv) f ai d{t, •) G L 6 (R) for all 5 > \/\d — 1|; in particular, f ad (t, •) G L X (R) n L 2 (R). 


4 














MvN—Kernel versus modified MvN—Kernel Role of the shift parameter for the modified MvN—Kernel 




Figure 2: Left: Comparison of the MvN kernel f^ vN (t , •) and the modified MvN kernel / a> (t, ■) for 

d — 0.25, a — 0.1. Right: Modified MvN kernel f ad ( t, •) with d — —0.25 for different values of the shift 
parameter a for f = 15. 


Proof, (i) and (ii) are obvious; (iii) is based on the fact that for all t G R, by a Taylor expansion. 


lim 

s — y —oo 


I /m 0' s ) 

IsK -1 


lim 

s —> — oo 


|(fl — s) d — (f + a — s) d 

IsK -1 


lim 

s—oo 


| dt(a — s) d 1 

Is |rf—i 


\dt\. 


(3.4) 


(iv) From (iii) follows that it suffices to check for which S > 0 the integral J 1 ! ds < oo for some 
N < 0. This is exactly the case for 5 > 1/ \d — 11. □ 

For any two-sided subordinator S = (Sf)^^ with E(S^) < oo, by substituting the MvN-kernel/^^f, •) 
in J2.31 with f a ,d{f> •) for d < 0 and a > 0, we define the fractional subordinator S a,d as 


S“' d = [ f a , d (t,u) dS u , tGR. (3.5) 

Jr 

Proposition 3.4. Let a > 0 and d < 0. Further denote by S — (St)teR a two-sided subordinator satisfying 
E(Sj) < oo. Then S a t ’ d as defined in ( 13.5b exists for all t G R as limit in probability and as limit in the L 2 (Cl)-sense. 

Remark 3.5. There is some freedom on the choice of an appropriate kernel. 

(a) For comparison, recall from above that the MvN kernel is for d G (—0.5,0) and f £ E defined as 


1 _ 1 fj PCO 

f^ vN (t,s) = r(1 + rf) (gd(t-s) -grf(-s)) = Y(l + d) TsJs 1 (0,t)( v )( v ~ s ) ddv ' t>seR, 

see ESI Lemma 1.1.3]. This kernel can generate fBm andfLm on 1R for symmetric driving processes. 

(b) The approach in Q to construct a FICARMA process driven by a subordinator ivould suggest to use the kernel 

-1 d 


_ 2 ^ rOO 

T(1 + d) ~ds Is t W)(v)™n(a d ,(v-s) d )dv / t,s e E, 


for d G (—0.5,0) and some a > 0. For us it was computational advantageous to bound g d and not the 
integrand in the above representation. 

(c) j[l4lf suggest a tempered fractional kernel of the form 

e -A(t-8)+( f _ s) rf _ e -K-')+(- s )d + d Vf t,s G E, 
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for d G (—0.5,0.5) and A > 0, which also has an integral representation using tempered fractional integrals, 
see j[l4l Definition 2.1] for details. It belongs to U for all p > 1. The autocovariance function of the 
increments of the corresponding tempered fractional process decreases for small lags algebraically, but for large 
lags exponentially. It is therefore called a semi-long range dependence model. 

(d) firm generalize fractional Levy processes by allowing the kernel to be regularly varying, which allows to prove 
functional central limit theorems for scaled Ornstein-Uhlenbeck processes driven by such generalized fractional 
Levy processes. 

3.1 Sample path properties 

The fractional subordinator S a,d has smooth sample paths in the following sense. 

Proposition 3.6. Let d G (—0.5,0) and a > 0. Further let S = (St)teR be a two-sided subordinator satisfying 
E(S 2 ) < oo. Then S n,d as defined in (13.5b has the following properties. 

d) {S a f d )t >o has a.s. continuous sample paths of finite variation. 

(n) {S a t ' d ) t >o constitutes a semimartingale (with respect to the filtration generated by S) with characteristics given 
by (S a ' d , 0,0). 

(iii) IfS has strictly increasing sample paths, i.e. characteristic triple (jsr 0 , vs) such that either 7 > 0 or v( 0 , 00 ) = 
00 , then (S a t ' d ) t >0 has strictly increasing sample path. 

Proof, (i) Recall the representation of the two-sided subordinator St — —S_j_lr f <o} + Sf l{t>o} f° r inde¬ 
pendent identically distributed subordinators S'■ 1 ' and S* 2i with So = 0. Define (Ft ■— c {St : t > 0}, hence, 
F~q is the sigma algebra generated by S' 1 ' 1 . Then {Sf d ) t >q is adapted to ((Ft)t>o- Now we can write 

s a4 = s a,d + J‘( a d_( a + t _ u ^ dSu/ t > 0 . 

Since the kernel function f a ,d{L') an d almost all paths of S have finite variation, we can apply partial 
integration to such paths, which yields 

Sf d = S a 0 ' d + d J^(a + t — u) d - l S u du, t > 0. 

This implies that ( S a t ' d )t>o has a.s. continuous and finite variational sample paths. 

(ii) Together with (i) Ifl6l Section II, Theorem 7] implies that (Sf d )t> 0 is a semimartingale with the given 
characteristics. 

(iii) This is a consequence of Prop. 1 of 0. □ 

3.2 Cumulant generating function and moments 

Let S be a subordinator without drift, then its cumulant generating function is given by In ]E(e (9S| ) = 
ip(O) = /“(l — e 6z )vs(dz), where 1/5 is the Levy measure of S. Proposition 2.6 of |I7| gives the cumulant 
generating function of S a,d as 

lnE(e 0S '' rf ) = InE^ exp \qJ fa,d{t,u) dS u j) = J n f Q C 1 “ e e ^ d(f ' u)z )i/ s (dz)du. 

The fc-th cumulant of S d,d is given by 

ASF = 


6 



For the fc-th derivative in 0 we obtain, provided the corresponding Levy moment exists 

K k {S a t ' d )= f z k v s {dz) f f k 4 {t r u)du^K k {S 1 ) I' f k d (t,u)du. 

J 0 J IR J ]R 


From this we calculate the mean and variance as 


E(Sf' d ) = E(Si) [ fad(t' id)du, 
JR 

Var(S?' rf ) - Var(Si) [ft 4 (t,u)du. 

J K 


3.3 Properties of increments 

Next we ask whether f a j serves its purpose in the sense that the increments of S a,d exhibit an appropriate 
dependence structure. We start with an auxiliary result. 

Lemma 3.7. Let a > 0 and d G (—0.5,0). Then the modified MvN-kernel f a4 as in ( 13.3b satisfies for t > 0 
/ r /«,A u) du = C + a 2d t - ^-(f + a) d+1 + -^(f + a) 2d+1 + c{t)t 2d+1 

with 


c = "“ +1 (jtT-mTt) mi ‘(‘) = J_f'[0-y) d -(-n) 


-,2 


dy, 


(3.6) 


and 


r(d + i) l 

limc(t) - r ( 2 d + 2 ) sin( zr(d + 0 . 5 )) + 2d+ 1' 


(3.7) 


Proof. By substituting y {u — a)/t we obtain 


f 

Jo 


a d — (a + t — u) c 


[ fldi^n) du = [ \{a-u) d -{a + t- uf 

j]R J —oo L 

= JL‘ U \( 1 ~ }L 7 1 ) - (~~T ~) 1 + + 


du 


t 2d+ 1 [ (1 — y) d — {—y) d dy + f a 2d — 2a d [t + a — u] d + [t + a — u] 2d du 

J —oo L J J o 


= t 2d+1 c{t) + a 2d t + [t + a- u] d+l * 


1 r. l2rf + l f 


o 2d +1 


[t + a — u] 


T(2d + 2) sin( 7 r(d + 0.5))' 


Further note that for the normalization constant 1/T(d + 1) in 13.1b we obtain 

/_ [{l-y) d + -{-y) d + \ dy = T(d + l) J ^ (/,f DN (l,y)) dy - r ( lf + D 

Consequently, 

c{t ) =_ r ( rf + 1 ) _ 

W T(2d + 2) sin(7r(d + 0.5)) 

F(d + 1 ) 

“ r(2d + 2) sin(7T(rf + 0.5)) J-a/t 
Since f_ a/t (1 — y) A — {—y) d dy —> 0 as t —> oo, the assertion holds. 


/ 

J —a/1 

(i -y) d -{-y) d 

f 

J —a/1 

(1 -y) d -(-y) d 


□ 
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Proposition 3.8. Let a > 0 and d £ (—0.5,0). Further denote by S = (Sf)i e jR a two-sided subordinator satisfying 
E(S 2 ) < oo. Let S a,d be defined in d3.5b , then the increments of S a,d have the following properties. 

(i) S a,d has stationary increments. 

(ii) Let r > 0 be fixed and s + r < t such that t — s — hr for some h > 0. Then the two increments Sff r — S n t ' d 
and S a s ’ d — Ss ,d of length r have covariance 


7 r (h) Cov 


( ca,d nCi/d nU,d rci,d \ 

s+(/z+l)r _ ^ s+hr ' ^s+r — J 


which satisfies 


jr^h) ~ Var(S"' d ) \d\a d r 2 (hr + a) d 1 , hoo. (3.8) 

Proof, (i) For n £ IN, to < h < • • • < t n and a\,...,a n £ JR we use the Cramer-Wold device and calculate 

Yh a i(St'+h ~ St’ d +h ) = J2 a i [ (ga,d(ti-l+h - u) - ga,d(ti + h - u))dS u 

i =1 i =1 • /R 

" [ 

J2 a i / (gnAti -1 -V)~ gaA l i “ v)) dS v 
1=1 711 


Z — 1 

n 


i=i 

where we have used the stationary increments of S. 

(ii) We introduce the notation St St — E(Sf) and S”' d ,/ R / fl/ d(h u ) dS u . For t,s > 0 we calculate 

Cov(S“' d , Sg ,d ) =Cov(S“' d ,S"' d ) = E (S a t ' d S a s ' d ) = i (E((S“' d ) 2 ) + E((S“' d ) 2 ) - E((S?' d - S a /) 2 
= i (E((S“' d ) 2 ) + E((S s fl ' d ) 2 ) -E((S?l d ) 2 )) . 

In the last step we have used that the increments are stationary. Furthermore, 

E((S“ /d ) 2 ) - Var (Sf d ) = Var (S a ) / / 2 d (t, W ) du, 

J IR, 

such that by the linearity of the covariance operator, 

7r(*) = Cov(S^ + 1 )r ,S s \ d r ) -Cov(S s ^ (;;+ 1 )r ,S^) -Cov(S^S^) +Cov(S^ r ,S s ^) 

- 5 (E((% d +1)r ) 2 ] +E((S«' d _ 1)r ) 2 ] -2E((S7;/) 2 ]) 

= ivar (Si) [ fa id ((h + 1 )r, u) du + f f 2 d ((h - l)r, u) du -if fl 4 (hr, u) du 
Now, using Lemma [3.71 we obtain 


7 r (h) = -Var(S 1 ) 


2a d 
d + 1 


hr 


+ a) + r ) d+1 + ((hr + a)~ r ) d+1 - 2 (hr + a) d+1 ) 

+ 2rf ! +1 (((fcr + a) + r ) 2d+1 + ((hr + a) - r) 2d+1 - 2 (hr + a) 2d+1 ) 
+ c(hr + r) (hr + r ) 2d+1 + c(hr — r) (hr — r ) 2d+1 — 2 c(hr) (hr) 2d+1 
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where c(t ) is defined as in ( 13.6b and according to ( 13.7b converges for t —> oo to a positive constant, which 
we denote by c. Consequently, a Taylor expansion and the fact that lim/,_ >co c(hr ± r) / c(hr) —> 1 gives for 
/z —> oo 


*<*> - ^ Var <*> [- m (ftr+ fl)d+1 ( +1)d ra + 


(hr + a) 2 


O 


(hr + a) 2d+1 
2d +1 


(2d + l)2d 
r 2 

+ c ( (2d + l)2d ——rj + O 


(hr) 2 


(hr + a) 2 
1 


(hr + a) 4 
b O 


(hr + n) 4 


(hr) 1 


Var (Si) (—d)a d r 2 (hr + a) 


\d -1 


□ 

Remark 3.9. fflj T/ic result (13.81 implies that asymptotically the increments of S a,d are positively correlated. 

(b) The increments of S a,d do not have long memory in the standard sense, since the covariance function is integrable. 
However, it decreases algebraically, and for d close to zero we approximately obtain long memory. This is in analogy 
to the asymptotic rate of decay of the modified CARMA kernel g a/ j in the case of a subordinator-driven CARMA 
processes in Section 8 o/@. 


4 The fractionally integrated COGARCH process 


Taking the subordinator S — [L, L ' l) > and its fractional subordinator S a,d with the modified MvN kernel 
( 13.1b as driving process in 11.31 1, we can now define the fractionally integrated COGARCH(l, 1) process. 

Definition 4.1 (FICOGARCH(l, d, 1)). Let ocq, ay, fi > 0 and d £ (—0.5,0). Assume L to be a Levy process with 
E(L 4 ) < oo. Define the price process G with some initial value Gq and 

dGt — Vt-dLt t > 0, (4.1) 

where the squared volatility (cr 2 ) t >o is given as the solution of the SDE 

da} = —f>\(a} — a 0 ) dt + cc^a} dS a t ,d . (4.2) 


The process S a,d is the fractional subordinator defined in 13.51 1 with modified MvN kernel 13.31 1 and subordina¬ 
tor S [L,LY d \ where [L,Lf D '> is the discrete part of the quadratic variation of L. The model 14.1b with 
m is called fractionally integrated COGARCH (1,1) process with fractional integration parameter d or 
FICOGARCH(l,d, 1). The stochastic volatility model 14.2b is called FICOGARCH(l,d, 1) volatility process 
with fractional integration parameter d. 

We can state the solution of the SDE 14.2b explicitly. 

Proposition 4.2. Consider the FICOGARCH( 1, d, 1) volatility process as in 14.2b . Then for almost all sample paths 
the pathwise solution of the SDE 14.2b with initial value a} > 0 is given by 


with 



t > 0, 


X t = fit - 0LiS a t ,d , t > 0. 


(4.3) 
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Proof. Since by Proposition 13.61 (1) (S n,d ) t >o has a.s. continuous and finite variational sample paths, also 
(X) f > 0 has and we can apply integration by parts to obtain 

da} — ^a} + otofi J e Xs ds^j d(e~ Xt ) + e~ x, d ^Oq + J e Xs ds 
— a}d[ — fit + oiiS a t ' d ) + e~ Xt e Xt aof\dt 


= —fii (af — oc 0 )dt + ociafdS' 


a,d 


Thus, (14.31 satisfies the SDE (14.21 . 


□ 


In Figure [3] we compare the sample paths of two FICOGARCHfl, d, 1) processes with respect to dif¬ 
ferent choices of d £ (—0.5,0). The left column shows the price, return and volatility process of a 
FICOGARCH(l, —0.01,1) process along with the sample autocorrelation function (acf) of the volatility 
process. The right column depicts the same quantities of a FICOGARCH(l, —0.40,1) process, which is 
driven by the same Levy process. From Proposition 13.81 we would expect a slower decay of the acf of the 
volatility process for larger values of d. This is confirmed by the two sample acfs in the bottom row of 
Figure [3] For simulating the fractional Levy process we approximated the process by the corresponding 
Riemann sums as explained in Ifl2l Section 2.4]. 

4.1 Stationarity of the FICOGARCH(l,d,l) 

To formulate the stationary version of the FICOGARCH(l, d, 1) model, we have to extend it to the whole 
of R. Consequently, in Definition 14.11 we choose a two-sided Levy process (Tf)fgR as hi <12.2b . Then the 
stochastic differential equation (14.2b has solution 


at — e 


X ‘ (°o + a opi e Xs ds^, t £ R, 


(4.4) 


X f - fit — cc\S a t ,d , f £ R. 


Since the subordinator S = [L, LJ : , , ) has a modification with finite variational sample paths, also the 
fractional subordinator S a,d has one, and hence X. Moreover, by Proposition 13.81 (i), S a,d has stationary 
increments. Note that 

s a,d _ I /-oo (( fl + ( _s )) rf - ( a + - s)) d )rfS s + f { ° ((fl + (~s)) d - fl rf )dS s forf <0, 

| /°oo (( fl + (~ s )) d - ( a + ( l ~ s )) d )dSs + /J (a d - (a + (t - s)) d 'jdS s for t > 0. 

The following is an analog of Theorem 3.7 of (13]| for a fractional Levy process. 

Proposition 4.3. Let L be a two-sided Levy process with E(L^) < oo. Let d £ (—0.5,0) and a > 0. Then S a,d , 
defined in (13.5b with subordinator S = [L, L]* D ), has a modification that equals the improper Riemann integral 


Sf d = d f {(a + (~s)+) d ~ 1 ~ (a + (t - s)-^" 1 ) S s ds - da d ~ l J* S s ds, t £ R. 


Moreover, S a,d is continuous in t. 


(4.5) 


Proof First note that for every finite variational sample path of S a,d we can apply partial integration and 
obtain for t > 0 


S a / = lim 


(^{a(—s)) d — (a(t — s)) d ^jdS s J (fi d — {a + (t — s)) d ^dS s 

( a d - (a + t) d )S 0 - d [ (a + (t- s)) d-1 S s ds + ( a d - (a + t) d )S 0 
Jo 


10 



Figure 3: Simulation of the modified MvN-FICOGARCH price process G (top) with corresponding return 
process AG (second row), volatility process u (third row) and sample acf 7 ^ of a (bottom) driven by a 
compound Poisson process with rate 5 and normally distributed jump sizes with mean zero and variance 
one half. The model parameters are ocq — 0.0195, — 0.0105, f'>\ — 0.0513 and a — 1. The fractional 

parameter d is equal to —0.01 (left) and —0.4 (right). 
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(4.6) 


- lim ((a + (-u)) rf - ( a + (t ~ I 0) rf ') s i( 

u — y —oo \ / 

+ lim d/ ((a + (-s)) d_1 - (a + (1 - s)) d_1 )s s ds 
Now recall that by the SLLN ^4' E(L 2 ) as s —» — oo. Hence, the integrand in ( 14.6b satisfies for s —y — oo 

((« + (-s)) d_1 - (« + (f - s))^ 1 ) |S s | - E(L 2 )(d - lM-s)*" 1 -t- 0 . 

Thus the integral ( 14.611 converges finitely for u —> — oo. This yields 

S a t ' d — d J - s) d_1 — (a + (t — s)) d_1 ^ S s ds — d J (a + (f — s)) d_1 S s ds. 


For t < 0 we get 


^a,d _ 


d J (^(a — s) d 1 — (a + (t — s)) d 1 ^S s ds + d (a — s) d 1 S s ds, 

Continuity follows from dominated convergence. 


| £d,d ga,d | 

Proposition 4.4. Let S a,d be as above with d G (—0.5,0). Then lim —-p— = 0 a.s. 

V-t—co 1 1 — V\ 

Proof. Wlog assume that v < t < 0 for fixed t. Then 

I S a t ’ d - S a /\ <\d\ ( f (a + (v- s))^ 1 -(«+(*-s)) 

\J — OO 


d-1 


S,ds 


+ I (a + ( ~s)) d 1 - (a + (f - s)) d 1 S s ds 
Jv 

+ J ( fl +( —s )) rf 1 [h,L] s ds^ —: |d| (/i + /2 + /3)- 


To estimate J\ we set u t — s and use a first order Taylor expansion to obtain for |0| < 1 


Ji 


A/T /*°° 

< 77 - 1 / (fl + u + {v - f)) d_1 - (a + M) d_1 (u + |f|)di< 

f — V\ Jt—v 


\t — V I |f — 1 

A A P OO 

< 7--1 / (d - l)(i> - f)(fl + u + B(v - t)) d ~ 2 (m + |f|)du 

|t - v | df-p 

/ POO POO 

= M|d —l|0(|f| / (a + u + 9(v — t)) d ~ 2 du + / (« + u + 9{v — t)) d ~ 2 udn 

< M9(\t\(a -d(t-v) + t- v) d ~ x + (a-9(t-v) + t- v) d ~ 2 ^(t - u) 2 ) 

= Mfl(|f |(« + (1 - 0)(t - v)) d ~ l + ( a +(l-9)(t- v)) d - 2 ±(t - v) 2 ) 

—> 0, v —> — 00 . 

To estimate J 2 we set again u t — s and use a first order Taylor expansion to obtain for |0| < 1 


k 


< 


M 


\t — v\ 1 1 — v\ Jt- 
M\d — l||f| 


p CO 

/ (« + « — f) d_1 — (a + w) d—1 (u + |f|)dw 

Jt—v 


< 


\ t _ v \ — J f ^(a + u + 9\t\) d 2 (u + \t\)du 


□ 
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Finally, 


- + u + e\t\) d - 2 udu + \t\ j”{a + u + 0\t\) d ~ 2 di^ 

< (^^(t-v) d + \d-l\\t\(a + t-v + 0\t\) d -^ 

— t 0, v —> —oo. 


/3 < J (a+(—s)) d x s ds — J (a + u) d 1 (u+\t\)di 


< 


/ 0 M U d du + \t\ J^~ l \a + u ) d - 1 = ^(H - M) d+1 + + M - I/)* - A 


Hence, 




M — 


d + 1 


MAI 


For u — » —oo, the first term tends to 0 since d < 0. By (13.4b we have (a + |p| — |f|) d — a d ~ m\v\~\t\) d -\ 
which implies that J 3 also tends to 0 as v — t — 00 . □ 

Lemma 4.5. The integral f* e~^ x, ~ x ^ds exists a.s. for v —> — 00 . 

Proof. Note that for C G (0, p\ /ol \ ) we estimate, using Proposition 14.41 


f 

JV 


,-(x t -x, 


rt , s a,d — S a,d > 

'ds< J exp | — fi(t — s) — aj —^ —-j—|f - s||ds 


|f-s 

J v ex p{ “ (*~ s )(/i + ^ )}^s < ^ exp | - (f-s)(j 6 i -«iC)jds, 


|f-s 


and the last integral exists for c —> — 00 . 

Next we prove that the process (of )feR := (f<-opi flc e~( Xf_Xs )rf s ^ is stationary. 
Lemma 4.6. For all t\ < ■ ■ ■ < t m , m G IN, h G 1R, 

Proof. We use the Cramer-Wold device and calculate for fli,..., a ni G 1R and rn G IN, 


□ 


m rtj+h 


«o(8i 


m m /» 

EAU = E a i 

i=1 f=l J 


exp 


1 = 

m ft- 

J^ a ij ' exp j - /Si(f; - v) +et 1 (S a t ' d u - Sf/fj J-di> 


{ — Pi{U + h — s) + — S” ,d ) S-ds 


1=1 

m rt, 


J t t +h u v+h 


III pf f N 

e «»• y_ ro ex p {- h (u - v) +«i ( s a t f ~ s a /) j 


dz; 


i=l 

1 


E A' 


Ai £1 

by the stationary increments of S a,d . 

The main theorem follows now from the preceding results. 


□ 
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Theorem 4.7. Let the squared FICOGARCH(l,d, 1) volatility process (af) t& ^ be given as in (14.4b with erf = 
ocqPi e Xs ds independent of {L t )t> o, which satisfies E(L^) < oo, and parameters as in Definition 14.11 Then 
(erf) t>0 is strictly stationary. Moreover, the price process {Gt)t>o as defined in (14.H l has stationary increments. 

Our new approach applies immediately also to higher order COGARCH models. 

Remark 4.8. The COGARCH(p,q) process was introduced in 10. By using the same approach we can generalise 
the FICOGARCH(l, d, 1) model in a straightforward manner to its higher order analog. 

Let p and q be integers such that q > p > 1. Further let «o,«i,.. . ,a p G 1R,, ..., f>q G K, 0,Pq f 0 

and ftp+i — ■■■ — Diq — 0. Then we define the q x q matrix B and the vectors a and lq by 


( 


0 

0 


B = 


0 

\~Pi 


1 0 
0 1 


0 0 
~Pq -1 ~Pq -2 


-Pi 



/ M \ 




a 2 


0 

, a = 

tt-q— 1 

' h — 

0 


\ OLn / 


\1 


with B —f>\ if q — 1. Then for a two-sided Levy proceess (Lt)teR as dz 12.2b with E(L^) < oo, we define the 
squared z’olatility process (crf)t>o with parameters B, a,ao by 


erf = a.o + a T Y f/ t > 0, 


where the state process (Yf) t >o is the unique solution of the SDE 

d\ t — BYfdt + lq(o 0 + a T Y f )dS"' rf , t > 0, 

with inital value Yg, independent of L. Furthermore, S a,d is the fractional subordinator as defined in 13.5b . If the 
process (cf)f >0 is strictly stationary and non-negative almost surely, we define the FICOGARCH(p,d,q) process 
(Gf)f>o with parameters B, a, ocq and some initial value Gq as the solution of the SDE 


dGt — Ut-dLt, t > 0. 


The volatility process of the COGARCH(l,l) and also of the FICOGARCHf 1, d, \) is non-negative by definiton. This 
is not the case for the COGARCH(p,q) model. Therefore, conditions as formidated in 0 Theorem 5.1] have to be 
considered to assure non-negativity of the volatility process. 
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